AD-A14J  305 


UNCLASSIFIED 


PASSIVE  LINE-DF-SIOHT  STABILIZATION  FOR  AN  INFRARED 
SENSOR(U)  NAVAL  RESEARCH  LAB  HASHINQTON  DC  R  L  LUCRE 
28  JUN  84  NRL-S828 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS- 1963-A 


AD-A143  305 


NRL  Report  M2I 


Passive  Line-of-Sight  Stabilization 
for  an  Infrared  Sensor 

R.  L.  LUCRE 

Advanced  Concepts  Branch 
Optical  Sciences  Division 


June  28,  1984 


a_ 


NAVAL  RESEARCH  LABORATORY 
Washington,  D.C. 


Approved  for  public  release;  distribution  unlimited 


8'4  07  25  093 


I 


REPORT  DOCUMENTATION  PAGE 


ia  RtCORT  C  _  ASS'  F  CAT  lijN 

UNCLASSIFIED 


2t.  Dfc  Cl  ASS'  F  *CA  T  ON  DOWNGRADING  SCHEDULE 


4  PERFORMING  ORGANIZATION  REPORT  Nu^BERlS 

NRL  Report  8828 


6a  NAME  OF  PERFORMING  ORGANIZATION 

Naval  Research  Laboratory 


.  ADDRESS  <  it >  Staff  and  /IP  Code 

Washington,  DC  20375 


Ba  NAME  OP  FUNDING  SPONSORING  Ub  OFFICE  SYMBOL  19  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

ORGANIZATION  [  il(  applicable  *  j 


3  DISTRIBUTION  AvAilABil'Ty  of  report 

Approved  for  public  release;  distribution 
unlimited. 


5  MONITORING  ORGANIZATION  REPORT  NUMBER'S' 


7a  NAME  OF  MONITORING  ORGANIZATION 


7b  ADDRESS  'City  yiate 


Naval  Air  Systems  Command  AIR  340J 


8t  ADDRESS  '(  ity  ''late  and  /.IV  (  ode  1 

Washington.  DC  20361 


II  t  TlE  Innudf  Security  Ctauifn  atmn 

(See  Page  ii) 


12  PERSONAL  AUTHORiSl 

Lucke,  Robert  L. 


10  SOURCE  OF  FUNDING  NQS 


13a  type  OF  REPORT 

Final 


16  SUPPLEMENTARY  NOTATION 


1 3b  TIME  COVE  RE  O  14  DATE  OF  REPORT  ,Vr  Mo  Do> 

from  Sept.  82  to  Aug.  83  1984  June  28 


COSAT.  COOES 


18  SUBJECT  TERMS  (Continue  on  rece'te  it  neceuar y  and  identify  b\  block  number 

Stabilization 

Line-of-sight 


19  ABSTRACT  n  oiUinue  on  reiene  if  necenar\  and  identify  6\  b/orfc  numberi 

The  use  of  passive  gyroscopic  stabilization  in  the  vibrational  environment  of  a  P-3  aircraft  is  described 
theoretically.  The  conclusion  is  reached  that  even  a  fairly  heavy  instrument  (30  kg)  can  have  its  line  of 
sight  stabilized  to  a  precision  of  100  brad  or  better.,- 

5  Y . .  _> 

r  / 


i 


/ 

>  t  4 

yf  /L"  •> 


20  D*S  r  «•  Bi  t 

ON  availability 

ABSTRACT 

v  NC  t  ASS'*  *  ( 

"Nl  MiTED  X  SAME 

AS  RPT  O  T  ,C  USE  RS  .  .. 

22a  name  of  REV-  NG  Hie  'ND-viDuAL 

Robert  L.  Lucke 


DD  FORM  1473.  83  APR 


21  ABSTRACT  SECuRiT  y  CLASSIFICATION 


22b  teifphone  number  |22<  of f-ce  symbol 

In*  ludi  1  'l  l!  <  ■  Mfi 


(202)  767-3045 


E  O'  TON  OF  1  j  A  N  7  3  i  S  O  B  SO  i.  E  T  E 


Code  6520 


SE  CtjMH  v  (  i  ASSlE  I(.a  I  M.)\  i  ■(  THIS  PA.  A 


m 


i 


St  CUM  I  T  Y  CLASSIFICATION  OF  THIS  PAGE 


11.  TITLE  (Include  Security  Classification)  (Continued) 
Passive  Line-of-Sight  Stabilization  for  an  Infrared  Sensor 


St  CUR  IT  >  CLASSIFICATION  OF  THIS  PAGE 


CONTENTS 


INTRODUCTION .  1 

EQUATIONS  OF  MOTION  FOR  A  SPINNING  ROTOR .  2 

EQUATIONS  OF  MOTION  FOR  A  GIMBALED  ROTOR  .  3 

APPLICATION  OF  EQUATIONS  OF  MOTION  TO  A  REAL  SYSTEM  .  4 

General  Form  of  Driving  Term  and  Resonance  .  4 

Effect  of  Static  and  Dynamic  Imbalance  .  7 

Effect  of  Rotational  Vibrations .  8 

CONCLUSIONS  .  10 


/ 

M  J 


iii 


PASSIVE  LINE-OF-SIGHT  STABILIZATION 
FOR  AN  INFRARED  SENSOR 


INTRODUCTION 

A  need  exists  for  background  data  taken  under  realistic  conditions  to  support  the  development  of 
infrared  (IR)  surveillance  systems  for  military  applications.  We  propose  to  address  this  problem  by 
mounting  an  existing  IR  scanner  in  the  cabin  of  a  Naval  Research  Laboratory  (NRL)  P-3  patrol  plane. 
It  is  positioned  to  view  the  outside  world  through  one  of  the  P-3’s  forward  observation  windows  (the 
plastic  window  material  will  be  replaced  by  a  pane  of  zinc  selenide).  This  location  permits  the  camera 
to  be  manually  operated,  which  eliminates  the  need  for  a  complicated  and  expensive  remote  control 
system. 

As  with  any  camera  mounted  on  a  moving  platform,  the  instrument  must  be  stabilized  against  the 
P-3’s  vibration  environment.  The  sensor’s  31.8  kg  (70  lb)  weight  and  330  microradian  (/trad)  pixel 
size  require  the  stabilization  system  to  hold  a  heavy  package  on  a  fixed  line  of  sight  to  an  accuracy  of 
100  fi  rad  so  that  no  significant  resolution  degradation  occurs.  Problems  of  this  kind  are  usually  handled 
with  an  active  stabilization  system:  gyroscopes  attached  to  the  instrument  sense  changes  in  direction 
and  activate  feedback  loops,  and  torque  motors  keep  the  pointing  direction  constant.  Such  a  system  has 
the  advantage  of  being  fully  automatic,  with  no  operator  required.  However,  it  is  very  expensive,  espe¬ 
cially  for  three-axis  stabilization  to  an  accuracy  of  100  /trad.  The  main  reason  for  this  is  that  errors 
introduced  by  the  feedback  loop  and  torque  motors  (especially  the  latter)  require  tighter  tolerances  for 
the  mechanical  system  (particularly  bearings)  which  can  be  met  only  by  expensive,  high-precision 
machining. 

An  alternate  approach  is  to  use  passive  stabilization,  based  on  the  principle  that  a  gyroscope  will 
resist  rotation  about  an  axis  by  precessing  about  a  perpendicular  axis.  This  principle  is  familiar  from  toy 
gyroscopes  which,  when  supported  at  only  one  end,  resist  the  force  of  gravity  trying  to  rotate  them 
around  a  horizontal  axis  by  precessing  around  the  vertical  axis.  This  principle  is  applied  to  the  sensor 
stabilization  problem  by  attaching  a  gyroscope,  gimbaled  about  one  axis,  to  the  sensor,  so  that  it  is 
forced  to  move  with  the  sensor  in  one  direction,  but  is  free  to  precess  about  its  gimbal  in  a  perpendicu¬ 
lar  direction  (see  Fig.  1).  To  quantify  this  principle,  we  first  develop  the  equations  of  motion  of  a 
freely  spinning  rotor,  then  couple  the  rotor  to  a  sensor  package..  The  theory  of  gyroscopes  has  been 
exhaustively  studied  for  hundreds  of  years,  so  no  new  physics  need  be  developed;  the  discussion  given 
here  is  tailored  to  our  particular  problem.  The  small  angle  approximation  will  be  made  in  the  equations 
of  motion  where  appropriate:  0.1  mrad  pointing  stability  will  be  lost  long  before  the  small  angle  approx¬ 
imation  is  violated.  This  permits  a  simple  derivation  of  the  root-mean-square  response  angle  of  the  IR 
camera  to  a  given  sinusoidal  input  disturbance.  This  response  can  then  be  evaluated  for  the  transla¬ 
tional  and  rotational  vibration  environment  of  a  P-3.  We  find  that,  for  reasonable  values  of  system 
mass  and  gyroscope  rotor  angular  momentum,  the  stabilization  system  can  easily  counter  the  direct 
effects  of  vibrations.  This  leaves  only  mounting  system  flexure  and  bearing  imprecision  as  the  dom¬ 
inant  disturbers  of  steady  pointing. 


Manuscript  approved  March  16,  1984. 
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Fig.  1  —  Schematic  representation  of  gimbaled 
gyroscope  attached  to  sensor.  A  disturbing  force 
D  in  the  Y-Z  plane  results  in  precession  P  in  the 
X-Y  plane. 


EQUATIONS  OF  MOTION  FOR  A  SPINNING  ROTOR 

The  orientation  of  the  angular  nmmentum  vector  of  the  rotor  is  specified  by  the  polar  angles  9 
and  0;  0  is  the  rotor’s  rotation  angle:  0  -  angular  speed  (see  Fig.  2).  The  Lagrangian  (which  is  just 
the  kinetic  energy  in  this  case)  for  a  spinning  rotor1  is 

L  -  -y- (02  +  02  sin2  9)  +  -y  (0  +  0  cos  9 )2, 

where  /3  is  the  moment  of  interia  about  spin  axis  and  /t  (  -  /2)  is  the  moment  of  interia  about  the  two 
perpendicular  axes.  This  leads  to  the  equation  of  motion  for  9  : 

Ail  <\r 

-  I\9  +  (/3  -  / j)02  sin  9  cos  9  +  /30  0  sin  9  =  N9  =  external  torque. 

Similarly,  for  0  and  0  : 

~~ |/ 10  sin  2  9  +  /3(0  +  0  cos  9)  cos  #]  -  N+ 

-J~[/3(0  +  0  cos  0)J  - 


Fig.  2  —  The  direction  of  the  rotor's  an¬ 
gular  momentum.  L,  is  specified  by  the 
polar  angles  9,  <t>.  'P  is  the  rotor's  angu¬ 
lar  speed. 


*H.  Goldstein,  Classical  Mechanics,  Addison  Wesley  Publishing  Co  .  Inc  ,  Reading,  Massachusetts.  1965. 
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The  rotor  spins  on  nearly  frictionless  bearings,  so  A^~0,  and  is  gimbaled  so  it  can  turn  freely 
about  the  vertical  axis,  so  /V*=0.  Since  we  are  interested  in  the  stability  of  the  gyroscope,  we  need  to 
solve  the  equations  of  motion  only  for  smalt  angular  displacements  around  an  initial  position.  We 
therefore  choose  an  initial  position  of  9  —  ir/2;  this  means  sin  9  «  1,  cos  9  -  0,  and  9,<t>,  i/>  are  now 
independent  angles,  as  a  brief  inspection  of  the  Lagrangian  shows.  We  therefore  introduce  the  small- 
angle  variable  r)  : 

i)  —  0  —  sin  9  =  1  +  0(t)2),  cos  0  ~  — i j  4-  0(tj3). 

The  equations  of  motion  are 

lyV  +  (l)  ~  /i)<£2(-t>)  +  /3(/l0  -  Ng 

and 

^[/l0  +  Ijitjt  +  <^(-1?))(-Tj)J  -  0 

or,  to  first  order  in  17,  for  the  <t>  equation  we  have,  without  loss  of  generality 

ly<j>  +  /30  ( — tj)  ~  0. 


Note  that  /3 ifi  —  Lis  the  angular  momentum  of  spinning  rotor,  so 

</>  “  j- V. 

which,  when  substituted  into  the  tj  equation  gives,  to  1st  order  in  rj: 

Ip)  +  L  ■  y-ij  -  Nt 


or 


..  ,  L2  Ng 

'  Tp~T 


As  we  expected  from  the  small  angle  approximation,  t)  undergoes  simple  harmonic  motion  with 
natural  frequency  o»0  -  L/l\,  which  is  just  the  familiar  nutation  frequency  found  in  gyroscope  theory. 
Once  n  is  found  by  solving  this  equation,  <£,  and  therefore  also  <f>,  can  be  found  immediately. 


EQUATIONS  OF  MOTION  FOR  A  GIMBALED  ROTOR 

We  must  now  put  the  rotor  in  bearings  in  a  housing  and  mount  the  housing  (gimbaled  on  one 
axis)  on  the  package  to  be  stabilized  (this  is  shown  schematically  in  Fig.  1).  In  the  9(—  17  +  n/2)  direc¬ 
tion,  the  gyro  and  the  package  must  move  together,  so  the  small-angle  equation  of  motion  is 

ffn  +  —  jV„terna|,  (1) 

as  before,  with  l\  replaced  by  the  moment  of  inertia  of  the  whole  system  computed  about  the  stabiliza¬ 
tion  axis.  If  the  system  is  forced  to  move  in  the  17  direction,  the  gyro  will,  while  it  is  resisting  this 
motion,  precess  in  the  direction.  Due  to  the  gimbal  bearings,  this  precession  motion  is  not  coupled 
to  the  package  (as  the  9  motion  was),  so  the  equation  of  motion  is  the  same  as  before: 

IGj>  +  L{- 17)  -0^0-  -j-r\  (2) 

‘G 

once  again,  where  lG  is  the  moment  of  inertia  of  the  whole  gyro  (rotor  and  housing)  computed  about 
the  gimbal  axis. 


R.  L.  LUCRE 


For  the  case  of  a  constant  driving  torque  (DC  input),  we  have,  from  Eq.  (1),  since  y  =  0, 

1  .  K 

*  L‘ 

Stabilization  will  remain  good  as  long  as  <f>  <  1  radian  (larger  <t>  means  that  the  gyroscope  has  precessed 
until  it  is  nearly  parallel  to  the  stabilization  axis  and  can  no  longer  stabilize).  Now 

so,  if  L  «•  40  kg  m2/s  (which  we  expect),  a  DC  torque  of  1  N  ■  m  can  be  tolerated  for  about  40  s.  In 
general,  we  can  say  that  any  constant  torque,  applied  for  a  time  r  such  that  Nr  <  40,  will  not  exceed 
the  capabilities  of  the  stabilization  system.  This  is  important  because  the  electrical  cables  and  gas  line 
(for  the  cryostat)  attached  to  the  sensor  exert  an  approximately  constant  torque.  This  "attaching 
torque"  can  be  minimized  with  good  design,  but  should  still  be  expected  to  precess  the  gyro  through  a 
large  angle  over  time,  a  condition  that  will  require  operation  intervention. 

Combining  Eqs.  (1)  and  (2),  the  final  equation  of  motion  is  now 

...  L2  N*  m 

t}  +  -7-7-1?  -  -7-  •  (3) 

lTlG  lT 

APPLICATION  OF  EQUATIONS  OF  MOTION  TO  A  REAL  SYSTEM 
General  Form  of  Driving  Term  and  Resonance 

The  driving  term  on  the  right-hand  side  of  Eq.  (3)  represents  input  to  the  stabilization  system 
from  translational  and  rotational  vibrations  of  the  aircraft.  Coupling  of  rotational  vibrations  through  the 
bearings  on  which  the  system  gimbals  are  supported  is  negligible,  as  will  be  shown  later.  Coupling  of 
rotational  vibrations  through  viscous  dampers  mounted  on  the  gimbal  shafts  is  given  by 

|y)  *  y(fi  ~  v) 

I  1  /damping 

where  the  angle  p  describes  the  angular  motion  of  the  airplane,  so  p  and  p  are  the  vibrational  angular 
velocity  and  acceleration,  respectively.  The  damping  depends  only  on  the  difference  between  p  and  vj; 
if  these  two  angular  velocities  are  the  same,  the  viscous  dampers  have  no  effect.  Since  the  airplane’s 
vibrations  will  be  given  as  a  Fourier  spectrum  p  —  Poe""',  we  can  write  p  —  (1  / io»)p.  With  this  part  of 
the  driving  term  separated  out,  the  equation  of  motion  becomes 


i?  +  YV  +  T~T  * 
1 T'G 


4  e""'  +  y  e" 
I  i  w 


N/I  now  represents  that  part  of  the  driving  term  which  couples  translational  aircraft  vibrations  into  rota¬ 
tional  moments  of  the  package— this  will  be  calculated  shortly.  Writing  the  natural  frequency  of  the 
system  as  «o0  -  (L2/ lTfG)l/2,  we  now  have  the  generic  equation  for  a  damped  harmonic  oscillator  driven 
by  a  sinusoidal  input: 


V  +  y  1?  +  W(f)? 


The  solution  to  this  equation  is 


or  +  iwy 


4 
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Ordinarily,  we  are  interested  in  the  magnitude  of  17  (not  its  phase)  for  a  given  drive  term  D.  Hence  we 
write 


JL 

D 


1 

1 

'/2  =J_ 

2 

T 

2 

V 

T 

2 

O>0 

1  - 

c 0 

+ 

(0 

JL. 

£l>0 

wo 

0)Q 

R  . 


Quantity  R ,  the  response  normalized  to  unity  at  o>  -  0,  is  plotted  in  Fig.  3  for  the  case 
a»o  “  2rr  •  10  Hz  *■  63  rad/s,  y  -  1  s'1  (critical  damping  is  yc  —  2  a>o  =  126  s-1,  so  the  damping 
assumed  here  is  quite  small:  y  <  0.01  yc). 


Fig.  3  —  Normalized  response  of  passive  gyroscopic  stabilizer  with  resonant 
frequency  F„  -  10  Hz.  For  total  response  multiply  by  (2irF„)-J. 


As  we  shall  see,  the  driving  term  D  will  be  obtained  from  the  vibration  spectrum  of  a  P-3  and  will 
be  given  on  a  per  square  root  of  bandwidth  basis,  so  the  root  mean  square  driving  acceleration  over  a 

bandwidth  A/ is  (£>JA/)1/2,  and  the  total  over  all  frequencies  is  [j*0  D2df  J  *  .  Therefore,  the  total 
rms  response  of  the  system  is  given  by  T)rm,  -  [17  f] 1/2,  where 


We  anticipate  a  sharp  resonance,  y  «  tu0,  so  virtually  the  entire  contribution  to  the  integral  comes 
from  within  a  few  y  of  «o>  where  the  integrand  is  well  approximated  by 
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and  D(<o)  =  £>(«0)  is  a  constant  over  the  narrow  bandwidth  of  resonance.  Therefore,  using 

dot  —  lirdf 


Vr  = 


D2 

8  iratn 


r 

Jo 


d<» 


(<o  -  <o0)2  + 


D2 

4*>oy 


» 


I 


» 


Hence 


D(.<oq) 

2tD0Vy 


(5) 


This  analysis  assumes  D(w)  is  reasonably  constant  over  w,  and  may  not  be  valid  for  cases  where 
Diu)  »  D( co0)  for  some  <u.  This  occurs  in  the  P-3  for  propeller  harmonics  at  17,  34,  68,  and  136 
Hz,  which  superimpose  a  discrete  line  spectrum  on  the  vibrational  continuum.  The  contribution  from 
these  discrete  lines  will  be  calculated  separately  later,  and  found  to  be  negligible.  P-3  vibration  data  for 
the  forward  observer’s  station  are  shown  in  Table  1  (see  also  Fig.  4). 


» 


Table  1  —  Typical  P-3  Vibrations  for  the  Forward  Observer’s 
Station  (from  Lockheed).  Note  that  1  g  —  10  m/s2. 


Low-Frequency 

Continuum 

17  Hz 

34  Hz 

68  Hz 

Vertical 

Longitudinal 

Lateral 

Roll 

Pitch 

Yaw 

0.003  g/VHz 
0.001 

0002 

<0. 1  rad/sVHz 
<0.05 

0.013 

0.01  g 

0.002 

0.002 

<0. 1  rad/s2 
0.05 

0.02 

0.03  g 

0.002 

0.006 

1  rad/s2 
0.05 

0.02 

0.06  g 

0.01 

0.1 

10  rad/s2 
0.6 

0.35 

Fig.  4  —  Typical  vibration  spectrum  for  P-3,  forward  observer’s  station 
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Effect  of  Static  and  Dynamic  Imbalance 

The  dominant  driving  terms  are  due  to  imprecise  static  and  dynamic  balance  of  the  package, 
resulting  in  aircraft  linear  accelerations  coupling  into  rotational  motions  of  the  package.  For  static 
imbalance  (referring  to  Fig.  5a):  given  a  package  with  mass  M  and  radius  of  gyration  r  (so  /  *  Mr2), 
out  of  balance  by  distance  b,  the  torque  resulting  from  a  translational  vibration  of  the  support  bearings 
given  by  aelml  ( a  is  the  acceleration  per  VHz)  is 


D  -  N/I  - 


So,  from  Eq.  (5), 


2(t>0yfy 


b  «=  0.001  m, 
r  —  0.3  m, 
a  —  0.1  m/(sVHz), 

<u0  —  2ir  •  10  Hz  -  63rad/s,  and 
y  —  (1  s)-1 


V,ms~  10  firad 

which  is  well  within  the  100  /a rad  pointing  stability  required.  Note  that  a  «  0. 1  is  a  generous  estimate; 
it  is  a  factor  of  2  or  more  higher  than  actual  P-3  vibration  data. 

A  similar  situation  exists  for  dynamic  inbalance  (referring  to  Fig.  5b).  For  a  differential  linear 
acceleration  between  the  support  bearings  of  ae'u\  the  torque  on  the  package  is 

N  -  Lima  •  ~  —  cma 


=  cma  e'“'. 


so  the  driving  term  is 


D-  N/I  - 


cma  ca  m 
Mr2  r 2  M 


-g-0.02, 

c  —  0.5  m,  and 
a  —  0. 1  m/  (sVHzT 


7}  —  80  /a  rad 


► 
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which  is  as  much  as  can  be  tolerated.  This  shows  that  achieving  c  •  {ml M)  <  0.01  m  is  desirable  for 
dynamic  balance,  but  we  should  again  note  that  the  value  of  0.1  for  a  is  too  high,  especially  for  the 
differential  acceleration  of  the  bearings. 


M  =  MASS  OF  ENTIRE  PACKAGE 
(SENSOR  AND  GYRO) 


1 

^  ROTATION  AXIS 

\ 

SUPPORT 

BEARING 

CENTER  OF 
GRAVITY 

1  MISSES  C.G. BY 

DISTANCE  b 

(a)  static  imbalance 


I 
I 
I 
I 
I 

l 
I 
I 

H 

(b)  dynamic  imbalance 

Fig.  5  —  Geometry  of  (a)  static  and  (b)  dynamic  imbalance.  In  (a),  rota¬ 
tional  motion  of  the  package  is  induced  by  a  translational  acceleration  of 
the  support  bearings  perpendicular  to  the  page,  while  in  (b)  a  difference  in 
the  acceleration  of  the  two  bearings  is  required. 


Effect  of  Rotational  Vibrations 

Rotational  vibrations  of  the  airplane  are  coupled  to  the  package  through  the  bearings  on  which  it 
is  mounted.  According  to  bearing  manufacturers  (e.g.  Fafnir  Corporation),  the  torque  transmitted 
through  a  bearing  is  nearly  independent  of  load  and  angular  velocity  (unless  they  are  very  high),  and 
can  be  expressed  by 

N  -  N0  Sgn  (0  -  it) 

where  NQ  (~  a  few  hundredths  of  a  nt  •  m)  is  a  property  of  the  bearing,  and  Sgn,  the  sign  function,  is 
±1  depending  on  the  sign  of  its  argument.  If  good  stabilization  is  assumed,  i.e.,  rj  «  /3,  then  tj  can 
be  neglected  and  all  we  need  to  know  about  the  above  expression  is  how  long  the  torque  will  have  the 
same  sign.  An  order  of  magnitude  estimate  of  this  time  is 

T  ^rms/ firms' 
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i.e.,  the  rms  value  of  p  divided  by  the  rms  value  of  its  rate  of  change.  Since 

firms  “  [J/32  df^ 

1  • 

and,  from  fi  —  ~P, 


we  can  find  t,  once  the  angular  vibration  spectrum  of  the  aircraft  is  known. 


For  a  P-3,  the  continuum  spectrum  of  p  is  reasonably  constant  over  the  frequency  range  of 
interest,  so  it  factors  out  of  the  integrals;  t  can  be  easily  estimated  to  be  <  0.1  s,  using  the  low- 
frequency  cutoff  fc  ~  0.1  Hz  discussed  below.  Combining  this  with  the  low  transmitted  torque  of  the 
bearing,  N0  ~  0.01  nt  •  m,  shows  that,  since  N0t  ~  10-2  «  40,  this  torque,  caused  by  rotational 
vibration  of  the  airplane,  represents  a  completely  negligible  disturbance  to  the  stabilization  system. 


For  a  given  component  of  the  P-3’s  discrete  vibration  spectrum,  fi  —  fio  e'“'  and  an  adequate 
approximation  to  the  driving  torque  can  be  obtained  by  noting  that  Sgn  (p)  =  Re(e'"0,  so  that  the 
driving  term  becomes 
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The  response  to  this  drive  is  therefore 

D  __  N0/ 1 

^  tog  —  to2  +  itoy  totf  —  to2 

unless  to  ~  (op.  The  first  discrete  frequency  in  a  P-3  is  o>  -  2n  ■  17  Hz  —  107  rad/s.  As  long  as  o>0  is 
kept  somewhat  smaller  than  this,  say  <o0  =  2ir  •  10  Hz  -  63  rad/s.  Then,  using  N0  <  0.01  nt  •  m, 
/  =  10  kgm2,  we  find 

T)  <  1  /xrad. 

The  other  and  more  important  coupling  of  aircraft  rotational  vibrations  to  the  package  is  through 
the  rotational  dampers  used  to  limit  the  system  response  to  resonance.  Equation  (5)  suggests  that  y 
should  be  large  (heavy  damping)  to  hold  down  the  amplitude  of  the  resonant  peak,  but  Eq.  (4)  shows 
that  a  large  y  increases  the  rotational  vibration  driving  term,  which  is  especially  to  be  feared  at  low  fre¬ 
quencies.  A  balance  must  be  found.  Setting  D(to0 )  —  y fi/too,  we  find  from  Eq.  (5)  that  the  resonant 
response  to  continuum  angular  vibrations  is, 

_  _  _  Jy  fi . 

Itoo'Jy  2a>o 

Using  fi  <  0.1  rad/s2  from  P-3  vibration  data  and  a>o  =  60  rad/s,  y=l  s_l  gives 

17  =  15  /zrad. 


The  low-frequency  response  is 


V 


(yP/w)2 
—  to2)2  +  to2y 2 


df 
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where  fc  is  a  low-frequency  cutoff  below  which  the  continuum  vibration  spectrum  does  not  extend.  A 
conservative  choice  is  fc  “0.1  Hz.  This  choice  is  motivated  by  the  fact  that  a  sinusoidal  angular 
acceleration  of  P  -  0.1  rad/s2  would  require  an  angular  motion  of  p  “  p/to2  =  0.3  rad  ~  15®,  which 
does  not  happen  on  a  P-3.  Thus 
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which  shows  that  y  =  1  s  1  is  a  reasonable  compromise.  (Actually,  y  will  probably  be  a  few  times 
bigger  than  the  1%  of  critical  assumed  here,  but  fc  will  be  larger  than  0.1  Hz  in  compensation.) 

The  vibration  environment  in  the  P-3  is  composed  of  a  continuum  part,  which  is  fairly  constant 
over  the  range  from  0  to  500  Hz,  plus  a  discrete  spectrum  of  propeller  harmonics,  of  which  the  dom¬ 
inant  components  are  at  17,  34,  68,  and  136  Hz.  At  the  forward  observer’s  station  (depending  on  how 
the  propellers  are  synchronized)  vibrations  at  these  frequencies  may  exceed  the  continuum  by  a  factor 
of  10  at  17  and  34  Hz,  and  as  much  as  100  at  68  and  136  Hz.  A  typical  example  of  this  vibration  spec¬ 
trum  is  shown  in  Fig.  4. 

The  stabilization  system  will  be  shock  mounted  in  the  P-3.  Therefore,  the  vibration  spectrum 
input  to  the  stabilization  system  will  be  the  P-3  vibration  spectrum  multiplied  by  the  transmissibility  of 
the  shock  mounts,  which  also  (obviously)  function  as  vibration  isolators.  The  transmissibility  of  the 
shock  mount  system  is  given  in  Fig.  6,  which  shows  that  the  discrete  frequencies  are  suppressed  by  fac¬ 
tors  comparable  to  their  enhancement  above  the  continuum  and  hence  are  unimportant.  This  is  done 
at  the  price  of  multiplying  a  section  of  the  low-frequency  continuum  by  the  resonance  peak  of  the 
mounting  system— about  a  factor  of  5.  This  presents  a  problem  only  if  it  coincides  with  the  resonant 
frequency  of  the  stabilization  system.  It  also  shows  that  the  resonant  frequency  of  the  shock  mounting 
system  (about  5  Hz)  must  be  kept  reasonably  separate  from  the  resonant  frequency  of  the  stabilization 
system  (about  10  Hz). 


Fig.  6  —  Transmissibility  of  typical  shock  mounting  system,  F„  -  5  Hz 
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CONCLUSIONS 

The  foregoing  analysis  shows  that  with  proper  balance  of  the  package  and  care  taken  in  preparing 
the  shock  isolation  system,  the  passive  approach  to  stabilization  can  easily  handle  the  vibrational 
environment  of  a  P-3.  The  stabilization  system  itself  does  not  introduce  unacceptable  errors.  This 
leaves  only  the  mechanical  tolerances  of  the  mounting  system:  if  a  shaft  is  mounted  on  two  bail  bear¬ 
ings  50.8  cm  (20  in.)  apart  and  each  ball  bearing  has  25.4  gm  (0.001  in.)  of  imprecision  (this  is  typical 
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for  standard  bearings,  higher  precision  bearings  are  available  only  on  special  order),  the  shaft  could  be 
expected  to  have  an  extreme  wobble  of  ±50  jmd.  Combined  with  a  small  amount  of  flexure  in  the 
shaft  itself,  this  shows  that  the  main  constraints  on  pointing  stability  are  mechanical  in  nature  and  can 
be  solved  by  standard  design  techniques. 


